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Abstract 

We study interactions between the dark solitons of the parametrically driven nonlinear 
Schrodinger equation, Eq.([T]). When the driving strength, h, is below y^ 2 + 1/9, two well- 
separated Neel walls may repel or attract. They repel if their initial separation 2z(0) is larger 
than the distance 2z u between the constituents in the unstable stationary complex of two walls. 



They attract and annihilate if 2z(0) is smaller than 2z u . Two Neel walls with h lying between 



\/j 2 + 1/9 and a threshold driving strength h sn attract for 2z(0) < 2z u and evolve into a stable 
stationary bound state for 2z(0) > 2z u . Finally, the Neel walls with h greater than h sn attract and 

P. 

annihilate — irrespective of their initial separation. Two Bloch walls of opposite chiralities attract, 

•i-H . 

while Bloch walls of like chiralities repel — except near the critical driving strength, where the 
difference between the like-handed and oppositely-handed walls becomes negligible. In this limit, 
similarly-handed walls at large separations repel while those placed at shorter distances may start 
moving in the same direction or transmute into an oppositely-handed pair and attract. The colli- 
sion of two Bloch walls or two nondissipative Neel walls typically produces a quiescent or moving 
breather. 
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I. INTRODUCTION 



This paper deals with the parametrically driven, repulsive nonlinear Schrodinger (NLS) 
equation: 

id t ^ + ±d 2 x ^ + ^ -\^\ 2 ^ = h^* -i^. (1) 

Here h is the strength of the parametric driving and 7 is the damping coefficient. In the 
absence of damping, i.e. when 7 = 0, this equation has the same stationary Bloch- and Neel- 
wall solutions as the (2:l)-resonantly forced Ginzburg-Landau equation and the relativistic 
Montonen-Sarker-Trullinger-Bishop model [l[ [3, [], Q] . However, unlike the Ginzburg-Landau 
and the relativistic model, where only one of these solutions (the Bloch wall) is linearly stable 
in their region of coexistence, the NLS equation exhibits Bloch-Neel bistability. Further- 
more, both domain-wall solutions of the NLS (also known as dark solitons or kinks) can 

n 

move with constant velocity, and the moving walls also turn out to be linearly stable [5j. 
This multistability is a rare phenomenon and leads one to wonder about the outcome of 
the soliton-soliton interaction. The problem that is of ultimate interest to a physicist, is 
formulated as follows: Given an initial configuration of many different walls, will the asymp- 
totic solution, as t — > 00, consist of predominantly Bloch walls, Neel walls or possibly some 
other, more complicated, structures? Will the walls form widely spaced, equidistant lat- 
tices, or tightly bound clusters? In this paper we make the first step towards answering 
these questions. 

When 7^0, equation (CQ) still has a dark soliton solution (in the form of the Neel 
wall); this solution was shown to be stable for all 7 jsj. Our analysis of soliton interactions 
will naturally include the case of nonzero damping, especially in view of the wide range of 
applications of Eq. ([I]) . 

Equation (TjQ) was indeed derived in a broad variety of physical situations. In fluid dy- 
namics, the repulsive ( "defocusing" ) parametrically driven NLS describes the amplitude of 
the water surface in a vibrated channel with large width-to-depth ratio . (The case 

of the small width-to-depth ratio gives rise to the attractive NLS {(], B, lsj].) The same equa- 
tion arises as an amplitude equation for the upper cutoff mode in a chain of parametrically 
driven, damped nonlinear oscillators. In the optical context, it was derived for the doubly 
resonant x optical parametric oscillator in the limit of large second-harmonic detuning 
Next, stationary solutions of Eq.([T]) with 7 = minimise the Ginzburg-Landau free 
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energy of the anisotropic Xl^-model. Here F = J FdX, where 

T = \(d x M) 2 - (1 + h)M 2 + iM 4 + 2hM 2 + F , 

and MpT) = (0, M y , M z ) is the magnetisation vector whose nontrivial components serve as 
the real and imaginary parts of the complex field *&(X) in ([ID: ^ — M y + %M Z . This model 
appeared in studies of stationary domain walls in easy-axis ferromagnets near the Curie point 
Nonstationary magnetisation configurations were considered in the overdamped limit: 
ty t = —5F/5ty* [12]. The damped hamiltonian dynamics — —iSF/S^* — 7^ provides a 
sensible alternative; this is precisely our equation (JTJ). 

Finally we note that Eq.([T]) also arises in a completely different magnetic context - 
that of a weakly anisotropic easy-plane ferromagnet in a constant external magnetic field 
j^]. Here, the external field (chosen parallel to the z-axis) forces the magnetisation to be 
almost homogeneous, and Eq.(TjQ) describes small deviations e\l/ = M x +iM y from the uniform 
magnetisation M w (0, 0, 1). (Here e is a small parameter.) In the XY-model, on the other 
hand, Eq.flT]) governs the magnetisation vector itself. 

The literature devoted to the interactions of dark solitons and kinks is quite extensive 
(although perhaps not as vast as for bright solitons and pulses). Kinks are known to attract 
antikinks in rea 
and 4 -theory 



-valued, single-component Klein-Gordon equations, such as the sine-Gordon 



13] , with and without damping terms [14 L The same is true for kinks of 
the real Ginzburg-Landau equations u t = u xx — V'{u) [14]. These results admit simple 
interpretations; the kink-antikink pair converges in order to minimise its total energy in 
the former case, and to minimise the value of its Lyapunov functional in the latter. Next, 
the dark solitons of the (undriven) NLS equations are known to repel [l5, [l(|. This can 
also be understood as an attempt to minimise the energy of the pair, at the expense of the 
gradient of the phase. Proceeding to the domain-wall solutions of the parametrically driven 
Ginzburg-Landau equation 

d t m = \d 2 x m + m- \m\ 2 m-hm\ (2) 

the Bloch wall and antiwall of opposite chiralities were shown to attract, while those of like 
chiralities repel 

,Q- (Here h < |, the region of stability of Bloch walls.) The Neel wall 
and antiwall always repel in their stability region {h > |) [17]. Unlike the case without 
driving, the interpretation of these interactions is not straight-forward though. 



As for the driven NLS equation, the behaviour of its Bloch and Neel walls is even less 
predictable at the intuitive level, while the variety of possible interacting partners is much 
wider due to the multistability. We will see that the interaction pattern is quite complicated 
indeed. The character of interaction (repulsion vs attraction) will depend both on the driving 
strength and the interwall separation. We will also show that it is influenced by stationary 
complexes of walls, stable and unstable. 

This paper constitutes the first part of our project; here, we restrict ourselves to inter- 
actions of the walls of the same type, i.e. Neel-Neel and Bloch-Bloch interactions. The 
analysis of the nonsymmetric situations, i.e. Neel-Bloch interactions, requires a different 
mathematical formalism and will be presented separately. (See the following publication 

Ho 

The dynamics of the repelling solitons is, in a sense, trivial: if the walls are initially at 
rest, they will simply diverge to the infinities. Less obvious is what the collision of two 
attracting walls will result in. The study of the asymptotic (as t — > oo) attractors arising 
in the parametrically driven NLS constitutes the second objective of our work. We will 
show that if the dynamics are dissipative, then, depending on the strength of the driving, 
colliding walls either annihilate or form a stable stationary bound state. In contrast to this, 
undamped collisions will be found to always produce a breather, a spatially localised, tem- 
porally oscillating structure. Depending on the initial conditions, the breather propagates 
or remains motionless, and in either case is found to persist indefinitely. 

The outline of this paper is as follows. The two fundamental solutions of Eq.([T]), the 
Bloch and Neel wall, are introduced in the next section. For the analysis of the Bloch-Bloch 
and Neel-Neel interaction we use the variational method, under the assumption of well- 
separated walls. The method is detailed in section II III and the resulting finite-dimensional 
systems are analysed in sections HVl IVHt and IVIIII Section [TV] deals with the Neel walls; 
section IVHI is devoted to the Bloch walls of the opposite chirality, while the like-chirality 
walls are examined in section IVIIII The conclusions of the variational analysis have been 
verified in direct numerical simulations of the full partial differential equation. The numerical 
simulations allow to advance beyond the limit of well-separated walls; in particular, we use 
this approach to examine the outcome of soliton collisions. We have allocated a separate 
section (section IVIj) to the Neel-Neel simulations while the Bloch-Bloch simulations are 
reported in sections IVHI and IVIIII along with the corresponding variational results. The 
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nontrivial attractors mentioned above — the stationary bound state in the case of dissipative 
dynamics, and the breather in the undamped system — are further investigated in sections 
IVl and IIXI respectively. Finally, the main results of this work are summarised in section |Xl 



II. DARK SOLITONS: PRELIMINARIES 

Dark solitons are localised patches of low intensity (\^\ 2 <C 1) in a high intensity back- 
ground (|\I/| 2 = 0{1)). The admissible backgrounds are described by the stationary, spatially 
homogeneous nonzero solutions of Eq.([T]), 

*jg t = a ± e«*, (3) 

where 



A ± = Jl±^/W^ 



r 



29 + = arcsin — , 29- = tt — arcsin — . (4) 
h h 

As one can easily check, for h < 7 the zero solution is the only stable background and so 
dark solitons can only exist for h > 7; this is the condition we are implicitly assuming in 
this paper. It can be shown that \ is stable for all values of h > 7 while , is always 
unstable. We will only consider dark solitons existing over the stable background so that all 
solutions obey |\1/| 2 — > A 2 + as \X\ — > 00. Accordingly, A stands for A + and 9 for 9 + for the 
rest of the paper. 

It will be convenient to transform Eq.fll]) so that the asymptotic solution is independent 
of h and 7. To this end, we set 

(X, t) = iAe i0 ij;(x, t), x = AX. (5) 

(We have also rescaled the spatial dependence for later convenience.) Under this transfor- 
mation, Eq.fJTJ becomes 

iifh + f ifrxx ~ A 2 |V|V + ^ + (A 2 - + - r) = 0. (6) 

This is the form of the parametrically driven, damped NLS that we will be working with in 
this paper. The stable background solutions of Eq.([6]) are simply V'flat = 
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Solutions to Eq.flHJ) with |i 
or nontopological [?/>(— oo) = 
is the Neel, or Ising, wall Q, Q, Q, 



1 as |x| — > oo can either be topological [if)(— oo) = — ^>(oo)] 
o o)]. Eg.® has two explicit topological solitons. The first 
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20|: 



iPn(%) = tanh(x), 



(7) 



named after the Neel wall in magnetism, which is a domain wall with the magnitude of the 
magnetisation vector vanishing at its centre. The second topological solution, which exists 
only for 7 = 0, is usually referred to as the Bloch wall [l, [3,1^, 4|: 



^b{x) = tanh(_Bx) ± iK sech(Bx). 



(8) 



(In magnetism, a Bloch wall is a domain wall connecting the two domains smoothly, with 
the magnetisation vector remaining nonzero everywhere.) In Eq.([8l), 



and 



B 



2vM 2 - 1 
A 



V4 - 3A 2 
A 



Ah 
l + h 



1 - 3h 
l + h' 



(9) 



The solutions obtained by multiplying ip^ and ipB by (—1) will be called antiwalls, or 
antikinks. 

The Bloch wall ([S]) exists in two chiralities. In the Bloch wall with positive imaginary 
part, the phase of the complex field ip(x) decreases (that is, the point on the unit circle 
moves clockwise) as x varies from —00 to +00. By analogy with the right-hand rule of 
circular motion, we refer to this wall as the right-handed Bloch wall. The phase of the Bloch 
wall with negative imaginary part increases as x increases (that is, the phase vector rotates 
counter-clockwise). This corresponds to a left-handed sense of rotation, so that this wall 
will be called the left-handed Bloch wall. The antiwall obtained by multiplying ipB by — 1 
has obviously the same chirality as its parent wall, ipB- Regardless of their chirality, Bloch 
walls only exist for h <\. 

In Ref.|5(, it was proved that the Neel wall is stable for all h > 7. The stability of the 
Bloch wall in its entire domain of existence was demonstrated numerically Q]. 

Examples of nontopological dark solitons are given by stationary complexes, or bound 
states, of domain walls. These can be formed by two dissipative Neel walls, or by a Bloch 
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and Neel wall — in the undamped situation. (See {5], Q.) In condensed matter physics, 
;hese nontopological solitons describe bubbles of one thermodynamic phase in another one 



221 ] . Accordingly, we will occassionally be referring to bound states of domain walls as 



solitonic bubbles. 



III. INTERACTIONS BETWEEN THE WALLS: THE METHOD 

In this paper, we will be paying special attention to the simplest situation where the 
interacting walls are initially at rest. Choosing the origin of the coordinate axis midway 
between two Bloch or two Neel walls, one can verify that the subsequent evolution will 
preserve this symmetric arrangement (see below). A symmetric pair of well separated walls 
can be approximated by a product function 

ip = ip 1 (x,t)ip 2 (x,t), (10a) 

where <fi and <p 2 represent the individual walls in the pair: 

ipi(x, t) = tanh [B(x + z)] — iKisech[B(x + z)], (10b) 
<f2(x, t) = tanh [B(x — z)\ + i/t2sech[£>(x — z)\. (10c) 

In Eq.f llOp . K\ and «2 are time-dependent parameters accounting for the deformation of the 
walls due to their interaction. The variable z = z(t) gives half the distance between the walls; 
without loss of generality, we take z to be positive. Finally, the constant B characterises the 
width of the walls (B = B for Bloch walls and B = 1 for Neel walls). 

Our analysis is based on the variational method. For this method, we substitute the 
Ansatz (flOl into the action integral that gives rise to Eq.flS]): 

S = J Ce 2 ^ dt, (11) 

where 



c = Re J - f i^,i 2 - f h 4 + h 2 + ^ w + m 2 ] + l ^W- m 2 \ - y 

(12) 



S 



In what follows we introduce a small parameter e 



~2Bz 



. Integrating off the explicit 



x-dependence in (j!2p . produces a finite-dimensional lagrangian 



C = T — V, 



(13a) 



where 



and 



T = -7ri(/ci + ft 2 ) [1 - 2(1 - Ki«2)(e - 2e ; 



2tt 



-^"( e - 4e2 ) [*i(l-«2) + *a(l-«i)] 



(13b) 



T ^ A 2 i3 2 + 8A 2 -12. 2 2 . 2A 2 . 4 4 , ( 



3B 



3B 



(AA 2 - A 2 B 2 - A)z - ^{A 2 - 1 



B 



+ A e 2 [ 2 ^ 2 (4B 2 - 7) + (30A 2 - 5A 2 B 2 - 12) (k\ + « 



+2A 2 (£ 2 - 11K*: 2 - 4A 2 (^ + «*) + 4A 2 (k 2 



)/%K 2 ] 



47T7 



+16e 2 z [2A 2 (1 - £ 2 ) + (A 2 B 2 - 6A 2 + A)(k\ + K 2 ) + 2A 2 K 2 K 2 ] +_i e (l-4e)(« 1 +K 2 )(l+K 1 K 2 ). 



(13c) 



[In ( 113bl) . the overdot indicates differentiation with respect to £.] The stationary action 
principle 8S = yields then the Euler-Lagrange equations for k 12 and 2. In deriving (113b|) - 
(I13cp we neglected all powers of e higher than e 2 . One can readily check that such terms 
would produce only higher-order corrections in the resulting equations of motion. 

Some comments must be made regarding the Ansatz ( fTOl) . After a variational Ansatz 
has been substituted in the corresponding field Lagrangian and the x-dependence inte- 
grated away, all dynamical variables (or possibly their combinations) can be grouped into 
canonically-conjugate pairs. (See 23j for review and references.) In particular, if the wave- 
form consists of two solitary waves and the separation between two constituents, 2z, is 
chosen as one of the variables, the conjugate momentum is the phase gradient. This fact is 



we 
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1 known in the case where the constituents are bell-shaped ("bright") solitons (see e.g. 



251]). That it remains true in a more general situation can be seen from an analogy 



with quantum mechanics where the momentum of the system is given by an eigenvalue of 
the operator \d x . (Since the eigenvalue has to be real, the operator acts just on the phase 
of the eigenfunction while the modulus is taken to be constant.) Another simple analogy 



arises if we write ip — \fP e%x ] this polar decomposition casts Eq.([T]) with h = 7 = in the 
form of equations of gas dynamics, where p is the density of the gas and d x x is its velocity, 
proportional to the momentum. 

In the case of bright solitons, each constituent soliton in the Ansatz is usually multiplied 
by e lk i x±z ) (see e.g. 24], 25]); then, the momentum conjugate to z is k. However, in the case 
of kinks this simple recipe would violate the boundary conditions at infinity. For this reason 
we have introduced the phase gradient by allowing the imaginary parts «i 2 sech [B(x ± z)] 
of the kinks to be variable. The imaginary parts decay as \x\ — » 00 and therefore making 
them variable is compatible with the boundary conditions. Below, we will show that k% and 
H2 are indeed momenta canonically conjugate to z. 

Another reason for introducing the phase gradient in this way is a simple interpretation 
of the finite-dimensional momenta k± and k 2 in terms of the original fields ip(x,t). Indeed, 
these variables coincide, modulo a numerical coefficient, with the field momenta of the two 
walls. The field momentum integral has the form 

P = \ J - W ) dx. (14) 

For the undamped NLS [equation with 7 = 0], the integral (fT4"j) is conserved. The field 
momentum of the Neel wall (171) equals zero, while that of the Bloch wall (jHJ) is given by 
Pb = T^Ko- [The same sign convention is used here as for Eq.Q — the top sign applies 
for the right-handed wall, while the bottom sign corresponds to the left-handed wall.] It 
is not difficult to check that the field momenta of the perturbed walls tpi and <y? 2 in the 
configuration fflOl) are -P[^i] = ttki and P[<£> 2 ] = —^^2- 

In symmetric situations that we are concerned with, the variables K\ and k 2 (both of 
which are conjugate to the same coordinate variable z) will eventually turn out to be related 
due to the equations of motion. The calculations simplify, however, if this relationship is 
established yet at the level of the lagrangian (fl3|) . that is, before the equations of motion 
have been derived. 

The relationship between K\ and k 2 depends on the type of walls being considered. If the 
initial configuration consists of two well-separated Neel walls at rest, or two quiescent distant 
Bloch walls of opposite chirality, we have a symmetry ip(x,0) = tp(—x,0). Since equation 
(J6j) is parity- invariant, the subsequent evolution will preserve this symmetry. Substituting 
the Ansatz fflO|) into ip(x,t) = ip(—x,t), we get K\ = k 2 - The initial value of n(t) should be 
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chosen close to zero for Neel walls and close to ±kq for Bloch walls [with the sign of k(0) 
depending on whether the wall placed on the right is right- or left-handed.] 

If the initial configuration consists of two quiescent Bloch walls of like chirality, the 
relationship between K\ and k 2 is not so trivial to establish. The reason is that equation (jH]) 
does not preserve the symmetry ip*(x, 0) = x, 0) of the initial condition. In this case we 
have to appeal to the field-momentum considerations rather than symmetries. Substituting 
the Ansatz flTD]) into the integral (fUj) . we find that P[ip] = 7r(/«i - K 2 )[l + 0(e" 2Bz )}. On the 
other hand, the total momentum of two like-chirality walls which are initially at rest and 
far away from each other, should be near =p27TKo- Therefore, for moving walls the quantity 
Ki — k>2 must not be different from =F2/t by more than 0(e~ 2Bz ). This is accomplished by 
letting Ki = =fko + q{t) and k 2 = ±k + q{t). When z is large, the perturbation q(t) is small 
— but not necessarily as small as e. 

The Euler-Lagrange equations corresponding to Eqs. (fTT!) . (fr3|) are 

dT_dV_ ±dT _ 2 dT = Q 
dn 8k dt dk dk 

and 

dT _ dV _ IdT _ ^dT = Q 

dz dz dt dz dz 

In Eq.f ll5al) . we assumed the situation of K\ — k 2 = k (that is, two Neel or two oppositely- 
handed Bloch walls). In the other case, i.e. when we have two Bloch walls of like chirality, 
n should be replaced by q in Eq.f ll5ap . The integrals T and V can both be found explicitly 
and simplified assuming wide separation of the walls. 

Finally, we note that the variational method is not well suited for the analysis of the 
Bloch-Neel interaction. The reason for this is that the Bloch and Neel walls have different 
widths, leading to terms in the Lagrangian which are not periodic along the imaginary axis 
on the plane of complex result, the Lagrangian cannot be obtained in closed form. 

(When the walls have the same width, the integrals can all be evaluated by integration along 
the rectangular contour with one side on the real axis and the other on the line Imi = vr/i3.) 
Consequently, we had to resort to a different approach for the analysis of the Neel-Bloch 
interaction (see paper II of this project 19].) The resulting phenomenology is also very 



different 19] . 
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IV. TWO NEEL WALLS 



In this case, we let B = 1 and K\ = K2 = k- In equations ( 1151) we discard products of 
powers of small quantities e and k for which there are larger counterparts. For instance, 
we drop a term proportional to e 2 n from an equation which already has a term en with 
a nonsmall coefficient. Here we keep in mind that in some cases the damping 7 and the 
difference \A 2 — || can be small parameters, hence we cannot drop e 2 K in favour of the term 
•yen or (A 2 — |)eK. With these simplifications, equations (I15ap - (115bj) take the form 



Kf 12 U 2 -^ 



K€ 



16 

+327T7eV + —A 2 K 3 e - Q4(A 2 - l)ne 2 , (16a) 
3 

7t(k + 27K) = 32A 2 e 2 
-24 ^A 2 - /t 2 e^ + 32(A 2 - l)e/t 2 . (16b) 

The subsequent analysis depends on the relation between the damping and driving in the 
system. 

A. 7 = 0(1) 

First we consider the case where 7 is not small. Assume, in addition, that A 2 is not close 
to |. (Later in this subsection we will explore the situation where \A 2 — || is small.) In this 
case the terms in the second lines of ( 11 6 aft and (116bl) are smaller than those in the first lines 
and can be neglected. The equations simplify to 

= 12 ^A 2 - |j «e, (17a) 
tt(k + 2 7 k) = 32A 2 e 2 . (17b) 

Since k is small, the variable e(t) varies very slowly, e/e ~ k. On the other hand, the variable 
«(t) will initially change on a much faster scale. Within the time At ~ 7 -1 it will "zap" 
onto the nullcline 

27T7K = 32A 2 e 2 , (18) 

after which the point (e, k) will be slowly moving along this parabola. According to Eq. (ll7al) . 
it will move towards greater e (i.e. the separation 2z will decrease) if A 2 > | and to smaller 
e if A 2 < |. That is, the walls will attract if A 2 > | and repel if A 2 < |. 

12 




20 40 60 80 100 
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FIG. 1: Vector field (a) A 2 < § (here A 2 = 1.2); (b) A 2 > § (here ^ 2 = 1.8). (c) gives a 

blow-up of a small neighbourhood of the origin in (b). 

Note that this criterion is consistent with results for the undamped, undriven case (7 = 0, 
A 2 = 1) where the dark solitons are known to repel jlol 

Since the above criterion has been derived under the assumption that A 2 is not very close 
to ~, the value A 2 = g provides only a rough watershed between the two types of interaction. 
In order to establish a more accurate borderline, we zoom in on a narrow strip along the 
line A 2 = | on the (7, v4 2 )-plane. Assuming that \A 2 — || is small (while 7 is not), Eq. fU7ap 
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should be replaced by 



, „ 9 4\ 64 o 64 9 
Tie = 12(A 2 --\Ke + —K 3 e-—K6 2 . (19) 



Substituting ([18]) into (|T9|) and dropping an 0(e 7 ) term, we obtain a one-dimensional dy- 
namical system 

256 / ~ 4 16 



- = — V^ 4 - 3 - y £ J £ • < 20 > 

When A 2 < |, e(t) tends to zero for all e(0). For A 2 > |, the system (1201) has a stable fixed 
point at 

€ = -(a 2 --). (21) 

16 V sy v ; 

This fixed point corresponds to a stable bound state of two Neel walls. 

When A 2 approaches |, the distance 2z = — lne between the walls in this bound state 
tends to infinity and so the bound state does not exist for A 2 smaller than |. This means 
that A 2 = | is, in fact, an accurate borderline between the two types of behaviour. For 
A 2 < |, two walls repel whereas for A 2 > |, the walls attract — except for A 2 close to |, in 
which case they form a stable bound state. 

We should emphasise here that our present conclusions pertain only to distant walls. In 
order to extend our understanding of the Neel wall dynamics beyond the limit of very small 
e, one has to analyse the dynamical system ffl6l) without neglecting any powers of e and k 
in it. This will be done numerically in section HV Dl below. We will show, in particular, that 
in addition to the stable bound state of two walls, there is also an unstable complex — at 
a shorter distance. 



B. \A l - || = 0(1); small 7 

When 7 is small, the system ffl6|) becomes 

vre = 12 ^A 2 - ne, (22a) 

n(k + 2 7 k) = 32A 2 e 2 - 24 ^A 2 - n 2 ez. (22b) 

For either sign of A 2 — |, the origin is the only fixed point in the system. Noting that the 
/t-axis (e = 0) is an invariant manifold flowing into the fixed point and calculating the index 
of the vector field (122]) along the semicircle e 2 + k 2 = R 2 , e > 0, with R < 1, we conclude 
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FIG. 2: Vector field j22J with 7 = 0. (a) A 2 < §; (b) ,4 2 > 



that the fixed point is a stable node for A 2 < | and a saddle for A 2 > |. Hence, for A 2 < |, 
all trajectories tend to the origin as t — > 00 [Figflja)], i.e., the walls repel. 

For A 2 > |, all trajectories fly away to infinity [Fig{T](b)] implying the attraction of 
the walls. One may wonder whether some trajectories coming from the lower half of the 
phase plane, could flow into the origin (tangent to the vertical axis.) A simple analysis of 
trajectories in a small neighbourhood of the origin shows that such behaviour is, in fact, 
impossible. Trajectories which seem to be approaching the point e = k = close to the 
ft-axis, do turn away when they are a small distance away from the origin. This is illustrated 
by Figfjjc) which gives a blow-up of a small neighbourhood of the fixed point. 

In the next subsection, we will zoom in on a small neighbourhood of A 2 = | to draw a 
sharper boundary between domains of the opposite types of interaction — similar to the way 
we have done it for large 7. We should also mention here that, as in the previous subsection, 
our current conclusions pertain to walls with large separations. 

The undamped case, 7 = 0, is exceptional; in this case, the entire K-axis is a line of fixed 
points — they represent pairs of infinitely separated walls moving at constant velocities. 
For A 2 < |, the points with k > are stable and those with k < unstable. Trajectories 
have the form of arcs starting on the lower semiaxis and flowing into points on the upper 
semiaxis: e(t) — ► 0, K,(t) — » > as t — » 00 [Figj2](a)]. The separation between the walls 
corresponding to the points on the upper semiaxis grows linearly: z — > — - (A 2 — |) K*t; 
that is, the evolution produces two Neel walls moving away from each other at a constant 
speed. 

In the undamped case with A 2 > |, the stable points are those with k < and unstable 
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those with k > 0. Since k > on the horizontal axis (where k = 0), trajectories starting 
on the upper semiaxis cannot cross to the negative-^ half-plane and have to escape to the 
positive-K infinity: k — > oo, e — > oo as £ — > +oo. The same is true for trajectories that 
have arrived from the negative-^ half-plane. Invoking the reversibility of the system (|22p 
with 7 = 0, i.e. making use of the invariance under t — > —t, k — > — k, completes the phase 
portrait [Figj2](b)]. 

To interpret the portrait, we note that k is proportional to z, the velocity of the wall 
[see Eq. (122ap ]. According to FigfS^b), two walls with zero initial velocities [k(0) = 0] will 
attract [e(t) — > oo as t — > oo]. The same applies, naturally, to the walls whose initial 
velocities are directed towards each other [n(0) > 0], and to walls with small outward initial 
velocities [«(0) < 0]. However, walls with sufficiently large outward initial velocities will 
diverge [e(t) — > 0], with the separation growing linearly: z(t) — > — | (A 2 — |) as t — > oo. 
(Here < 0.) The dashed line in Figj2](b) is a separatrix between initial conditions giving 
rise to the above two scenarios. 

Finally, it is worth mentioning that when 7 = 0, equations (fl5l) with K\ = K2 = k and 
B = 1 form a hamiltonian system (with a nonstandard bracket): 

1 dH . 1 dH 

z = —ir> K = — IT' 23 

TIT OK TXT OZ 

where 

r(e, «) = 1 + 2e(l - 14e + + lOe^ 2 ) > 0, (24) 

and the Hamiltonian 

H(e, k) = (1 - 3/i)/t 2 - -A 2 /t 4 + 4e/t 2 [(l - 3/i)z + 4/i] 

3 

+ -e 2 [3A 2 + (12 - 25A 2 )k 2 + UA 2 k 4 - AA 2 k 6 ] 
3 

+ 16e 2 zK 2 (5A 2 — 4 — A 2 k 2 ). 

This observation endows k with a simple interpretation of the momentum canonically con- 
jugate to z — in agreement with qualitative arguments in section IIII1 
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C. Small | A 2 - §|; small 7 



Finally, let us assume that both \A 2 — || and 7 are small. Here, equations (1T61) are to be 
replaced by 

( 4\ 64 o 64 „ . , 

vre = 12f A 2 - - Ue + y/t 3 e- y/te 2 , (25a) 

Tx(k + 2-iK) = ^t 2 + ^tn 2 . (25b) 

Note that we have neglected the term proportional to z = — | lne in the second line of (116bl) . 
When A 2 — | < 0, this is justifiable as this term is of the same sign as the last term in the 
right-hand side of ( Il6bl) and so it cannot alter the dynamics qulitatively. The validity of 
this approximation for A 2 — | > will be established below. 

The dynamics is influenced by nontrivial fixed points which arise as intersections of the 
nullclines of the system (1231) : 

/ 54 64 \ 

ne f 12/2+ —k 2 - —ej = 0, (26a) 

l^ e 2 + ^ £K 2 _ 27r7ft . = a (2 g b) 

Here we have introduced 

yU = A 2 -| (27) 

The nullcline ( I26b|) emanates out of the origin as a parabola: k — > Jp^e 2 as e — > 0; 
turns back and escapes to infinity asymptotic to the positive K-axis: k — > ^-e -1 as e — >• 0. 
Importantly, it has no points in the e > 0, k < quadrant. The nullcline (I26ap consists 
(apart from the vertical and horizontal axes) of a parabola lying on a side: e = |/i + |k 2 . 
Intersections of these nullclines are easy to determine and visualise, and the stability of the 
fixed points can be classified by straightfoward index arguments. 

When fi < 0, there is only one nullcline intersection in the positive-e halfplane. It is 
not difficult to check that this fixed point is a saddle. The origin is also a fixed point - 
a stable node. The portrait is in Figj3](a). Distant walls [more precisely, walls with initial 
conditions lying below the stable manifold of the saddle in FigJS^a)] repel: e(t) — > as 
t — ► 00. On the other hand, nearby walls (more precisely, those above the stable manifold) 
attract: e(t) — > 00. 

As /i grows through zero, another fixed point appears from the negative-e half-plane. An 
exchange of stabilities occurs at \i = 0: the nontrivial point appearing for \i > is a stable 
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node whereas the origin becomes a saddle. The nearby walls attract while distant walls form 
a stable bound state [FigJ3](b)]. As \x is increased further, the two nontrivial fixed points 
merge in a saddle-node bifurcation. When // > /j sra , the origin persists as the only fixed 
point in the system (an unstable one). All walls attract [FigJ3J(c)]. 

Note that the logarithmic term in (116bl) that was omitted in (125bp . becomes of the same 
order as the last term in the right-hand side of (116bl) only if zfi ~ 1, i.e. when e is as small as 
6 -2/m However, since the nullcline (126aj) is bounded from the vertical axis by the inequality 
e > YqI^j there can be no fixed points with e ~ e~ 2 ^. Therefore the logarithmic term can 
have no significant effect on the dynamics and its omission is totally legitimate. 

To find the threshold /i sn (7), we eliminate k between ( 12 6 aft and (126bft . The resulting 
equation can be written as 



7e — y^/x 37r 



7. (28) 



If /j is positive, the function F(e) in the left-hand side of (1251) grows to infinity as e — > yjr/i 



and e — > 00, and has a minimum at e = e m j n = y^g(31 + v457)/i. The minimum value is 
F(e min ) = F /i 3//2 , where F is a numerical coefficient. Consequently, Eq. (l28|) has no roots if 
H > (3tt/16F ) 2 / 3 7 2 / 3 « 0.30487 2 / 3 and two roots if < /i < 0.30487 2 / 3 . 

According to Eq. (l28l) . when 7 grows, the "smaller" (stable) fixed point tends to e = j^fi; 
this reproduces Eq.( l2Tl) of subsection IIV Al 

These considerations translate into the following criterion valid for small 7 and small 
differences A 2 — |. [We are also assuming that the imaginary parts of the perturbation 
of the walls are small initially, i.e. k(0) is small.] The direction of forces between the 
walls depends on their position relative to two stationary bound states, a stable (with 
the interwall separation denoted 2z s ) and an unstable one (with the separation 2z u ). The 
unstable complex exists for all A 2 < A 2 n while the stable one only for | < A 2 < A 2 n , where 

A 2 sn = i + 0.3048 7 2/3 . (29) 

In their region of coexistence (| < A 2 < A 2 n ), the stable complex has a larger separation: 
2z s > 2z u . As A 2 — > |, the distance 2z s grows to infinity. 

When A 2 < |, two Neel walls repel if their separation distance is greater than 2z u . If 
the initial separation is smaller than 2z u , the walls attract. It is natural to expect that 
the attraction should result in the annihilation of the walls, but this diagnosis is beyond 
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the scope of the variational approximation and can only be established in direct numerical 
simulations of the full partial differential equation (jf)]) (see below). 

Next, if | < A 2 < A 2 n , the walls with separations larger than 2z u will move towards the 
stable equilibrium [from inside or outside depending on whether z(0) < z s or z(0) > z s ]. 
The walls at shorter distances, 2z(0) < 2z u , are attracted to each other; they converge and, 
apparently, annihilate. 

Finally, for A 2 > A 2 n , the walls are attracted to each other. No complexes, stable or 
unstable, can be formed for these large A. 

Like in the large- situation, the undamped case (7 = 0) is exceptional in the small- 
region. In this case, the whole K-axis is a line of fixed points. When fi > 0, the phase 
portrait [FigH](b)] is similar to Figj2]^b). The walls which are initially at rest [/t(0) = 0], 
attract (e — ► 00 as t — > 00). As for the moving walls with large separations, their initial 
velocities are determined by k(0): 

(30) 

If the walls are initially diverging sufficiently fast, they will continue to do so, with 

«(t) -> < 0; z(t) -> Ufi + t. 

The portrait for fi < is in FigJH(a). Restricting our analysis to small |/c(0)|, we observe 
that walls with large separations repel and small separations attract. This is true both for 
initially quiescent and slowly moving walls. 



i(0) 



«(0) 



7T 



6/i+^ 2 (0) 



D. Numerical study of the fixed points; all 7 

As we have seen in section IIV[ distant walls with A 2 > | are attracted to each other. 
Whether the attracting walls are going to form a stable bound state or collide and annihi- 
late, depends on two factors: (i) whether a stable and unstable bound states exist for the 
corresponding 7 and A 2 , and if yes, then (ii) how the initial interwall separation compares 
with 2z u , the interwall separation in the unstable complex. For small 7, the region where 
the two complexes exist was found to be | < A 2 < A 2 n , with A 2 n as in Eq. ff29|) . In order to 
classify the outcomes of the wall attraction, we need to demarcate the corresponding region 
for nonsmall 7 as well. 
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The one-dimensional dynamical system (1201) has only one, stable, fixed point. In order to 
describe the unstable complex (which is bound tighter than the stable one, i.e. has a greater 
e), it is sufficient to retain the next, 0(e 7 ), term in Eq. lfTTjj) when we substitute <HM for k. 
However, despite capturing both the stable and unstable bound states, the system (fTHl) - (fTU|) 
will not produce any reasonably-accurate description of the saddle-node bifurcation where 
the two fixed points merge and disappear. One reason for this is that the term 32n'ye 2 K 2 in 
Eq. (jl6ap which was dropped from Eq. (fT9j) . becomes larger than e 7 when we substitute f fT8|) 
for k. Another problem is that Eq. (118j) keeps only the leading term in the k(e) expansion; 
retaining the next, 0(e 5 ), term in k(e) produces another term larger than e 7 in Eq. (1201) . 

Thus, in order to obtain an accurate variational description of the saddle-node bifurcation, 
we need to return to the original, nonreduced, vector field ffT75j) . Denoting f(e, n) and g(e, k) 
the right-hand sides of (I16ap and (I16bj) . the two fixed points arise as intersections of the 



nullclines 

/(e,«) = 0, g{e,K) = 0. (31) 

The nullclines have a common tangent when 

fe{e,n) = g e (e,K) 

Using a standard numerical routine to solve the system (l3TT) -( [32l) for the vector of unknowns 
(e, k, A 2 ), one can find the saddle- node value A 2 sn for each 7. 

The resulting curve A 2 = ^4^(7) is shown in Fig J7] below (the dotted line), where it is 
compared to the data from the numerical continuation of the wall complexes as solutions 
to Eq.([6]). The significance of the curve is that it provides a borderline between the two 
dynamical scenarios: For A 2 > A 2 n , two Neel walls attract each other and annihilate; for 
I < A 2 < A 2 n , the walls attract and annihilate if the initial separation 2z(0) < 2z u and form 
a stable bound state if 2z(0) > 2z u . 

Note that the above numerical result is valid both for large and small 7. For small 7, our 
numerical curve A 2 sn (j) is reproduced by the explicit formula fl29l) . 

We conclude this section by mentioning that the unstable fixed point exists for all A 2 < | 
for which there are Neel walls, i.e. for 1 < A 2 < |. This is a result of the numerical study 
of the system (13T1) with 7 varying from to 1 in steps of 0.01, and A 2 varying from 1 to | 
in steps of 0.001. The reservation that we should make here, however, is that when 7 tends 
to and, simultaneously, A 2 tends to 1, the coordinate k of the saddle point approaches a 
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value of order 1. This contradicts the assumptions we made in the derivation of the system 
ffT6]) and so the fixed point in this parameter region cannot represent any bound states of 
the walls. (The value of the K-coordinate is reasonably small only when 7 is greater than 
0.5 or when \A 2 — || is smaller than 0.01.) Consequently, the variational analysis cannot 
provide a trustworthy description of the small- distance dynamics of the walls in the 7 ~ 0, 
A 2 « 1 region. The direct numerical simulations of the full NLS equation ([6]) seem to be 
the only reliable source of information here. 

V. BOUND STATE OF TWO DISSIPATIVE NEEL WALLS 

In the previous section we showed that two bound states of Neel walls, a stable and an 
unstable one, exist in certain parts of the (7, /i)-plane. Here, we will reobtain these solutions 
numerically, demarcate their regions of existence and compare these to the corresponding 
variational results. 

We will be using the term "solitonic bubble" as a synonym for the wall complex. Treating 
bound states as independent soliton-like entities is meant to emphasise the significance of 
these stationary solutions as possible attractors in the phase space; also, it should reflect 
their different, nontopological, nature. In fact, when the complex is tightly bound, it looks 
more like a single entity than a wall doublet; see Figj5] 

Our approach here is based on the numerical continuation (path-following) of the sta- 
tionary bubbles as solutions to the ordinary differential equation 

fii> xx - A 2 ivf^ + ii> + (a 2 - i)r + z 7 (^ - r) = 0. (33) 

The variational approximation (flOl) with parameters found in subsection IIVCI was used as 
an initial guess for the numerical solution at small 7 and h ~ |; this, in turn, served as a 
starting point for our continuation. We classified the stability of the resulting stationary 
complexes as solutions of the full partial differential equation ([6]), by linearising about the 
solution and examining the spectrum of the linearised operator numerically. 

To present results of the continuation graphically, we use the energy functional 

E-yfa + W-l W+^h l}*. (34, 

Although the quantity defined in Eq. (l34|) is not conserved for nonzero 7, the energy can 
be used as a bifurcation measure for stationary solutions. [The field momentum (1141) is not 
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suitable for this purpose as it satisfies P = — 27P and so all stationary solutions with 7 7^ 
have the same, zero, momentum.] 

A typical pair of bifurcation diagrams (corresponding to the separate continuations in 
h and 7) is presented in Figj6j [Figj6] corresponds to 7 = 0.35 (main frame) and h = 0.5 
(inset); the continuations for other fixed values of 7 and h produce curves of similar shape]. 
Each diagram consists of two branches; the entire branch with higher energy is found to be 
linearly stable. The solutions on the lower branch are found to be unstable; the instability 
is associated with a (single) positive real eigenvalue. Numerical simulations of the full 
time-dependent PDE flHJ) show that when perturbed, the unstable bubbles decay to the flat 
solution. 

As we approach the termination point of the top branch, the separation of the walls 
becomes infinite. The termination point satisfies h ~ a/7 2 + 1/9 (i.e. A 2 « |) — therefore, 
this solution reproduces the stable bound state detected by the variational method (sections 
IIV Al IIV CI) . The solution on the lower branch corresponds to the unstable complex. We 
note that this lower branch can be continued all the way to h = 7, in agreement with 
variational results in section HVDl 

The walls pull closer together as we move towards the turning point along the upper 
branch (in the main frame of Fig. [6]). As we continue away from the turning point along the 
bottom branch, the distance between the walls continues to decrease, reaches a minimum, 
and then starts increasing. This nonmonotonicity can be explained in terms of the fixed 
points of the dynamical system (|2"5l) . 

The existence and stability of the bound states is summarised in FigJTJ This figure 
charts the (7, /i)-plane into regions of different type of interaction between the walls. In the 
shaded region and above it, distant walls attract each other. The attraction may result in 
the creation of a stable bound state at some finite separation; this happens in the shaded 
domain. [The formation of the bound state is illustrated in Figj8](c) below.] Above the 
shaded corridor, two attracting Neel walls are expected to collide and annihilate. Finally, 
in the region between the shaded area and the bisector line h = 7, two Neel walls at large 
separation repel. 

In FigJTJ we a l so display the variational results for the saddle-node bifurcation curve 
h sn ( / -f). The dotted line in the main frame gives the result of solution of the system of 
three equations (|3Tl) -( l32l) while the dotted line in the inset was plotted using Eq. (l29l . The 
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variational results are seen to be in good agreement with the outcome of the numerical 
continuation. 



VI. NUMERICAL SIMULATIONS: NEEL WALLS 



Conclusions of the variational analysis were verified in direct numerical simulations of 



the NLS equation (j5J). We used a split-step pseudospectral method [26(, with N = 512 
Fourier modes on the interval (— L/2, L/2) = (—40, 40). The numerical scheme is stable for 
the timesteps At < tt~ 1 (L/N) 2 ; accordingly, we set At = 2.5 x 10~ 3 . The method imposes 
periodic boundary conditions. The initial condition had the form ([TO]) where K\ and «2 
were set equal to zero. The results are shown in FigsJHJ where we have "zoomed in" on the 
solitons neglecting the flanks of the simulation interval. 

To verify the diagram of FigtZl we examined three sequences of points on the (7, /i)-plane: 
one just above the shaded region, one just below it, and one inside. The sequence just above 
the saddle-node line consisted of ten points 7 = 0.0,0.1, ...,0.9 with h = 0.4 + 0.757. For 
each pair of parameter values, we examined the initial separations 2z(0) = 3.0,4.0 and 5.0. 
The results were in agreement with the diagram of FigJTJ all simulations gave rise to the 
convergence of the walls. For 7 7^ 0, this resulted in the annihilation of the walls [FigJS^a)] 
while in the dissipation-free case, the convergence of the walls ended in the formation of a 
nonpropagating breather [FigJH](b)]. 



The sequence just below the line h = a/7 2 + 1/9 consisted of 7 = 0.0, 0.1, ...,0.9 with 
h = a/7 2 + 0.09. In the 7 = case, the walls separated by 2z(0) = 4 were observed to repel 
whereas those placed at a smaller distance 2z(0) = 3, attracted and formed a nonpropagating 
breather. These results are in agreement with the variational analysis [see FigJU^a)]. When 
7 is nonzero, the initial separations 2z(0) = 3.0 and 4.0 gave rise to the repulsion of the 
walls. In order to observe the attraction, we had to reduce the initial distance to 2,2(0) = 1.4 
for 7 = 0.1; to 2z(0) = 0.8 for 7 = 0.2 and to 2z(0) = 0.1 for 7 = 0.4. These observations 
are consistent with the fact that z u , the separation in the unstable complex, becomes smaller 
as 7 grows. 

Finally, the sequence of points through the dashed region in FiglT] included 7 = 
0.1, 0.2, 0.8 with h = | + 0.7847. For these values of h, the variational method pre- 
dicts a tightly bound stable complex of two walls and indeed, in all eight cases the walls 
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with initial separations 2z(0) = 3.0,4.0 and 5.0 moved towards each other and formed a 
stable bound state at a certain smaller distance [FigJHtc)]. 

Thus the simulations of the Neel wall interactions are in agreement with our expectations 
for distant-wall dynamics summarised in FigJTl 



VII. TWO BLOCH WALLS OF OPPOSITE CHIRALITIES 

Proceeding to the analysis of Bloch walls, we remind the reader that the Bloch walls exist 
only for 7 = and h < |; these will be our assumptions in this section. Here, we consider 
the interaction of two Bloch walls of opposite chiralities (i.e. with total field momentum 
P = 0). The corresponding lagrangian arises by setting Ki = n 2 = k and B = B in Eq. ljTBl) . 
The resulting equations of motion can be cast in the form of a hamiltonian system 

K = ~ — -HT> (35) 



TIT 8k ' TIT dz 

where r(e, k) is given by Eq. (1241) . and the Hamiltonian 

^,*) = — (« 2 -^) 2 + i^« 2 

+ — e 2 [7 -9h- 2(9 + 5h)n 2 + (15 + llh)^ - 4AV] 
3B 

+ l6A 2 e 2 z(2K 2 - k% - k a ), (36) 

with e = e~ 2Bz and k,q as in (Q. Trajectories of the system (1351) are simply level curves of 
the function H(e, «). 

When the interwall separation is large, k should be close to its unperturbed value: k « k 
or k ~ — Kq, where k,q is, in general, 0(1), but becomes small when h is near |. In the 
derivation of the variational equations, we assumed that the walls are well-separated and each 
wall is close to its unperturbed form. Therefore, only regions around the fixed points (0, k ) 
and (0, — k ) can be interpreted within the PDE (jSJ) with full certainty. (Trajectories outside 
those regions may also have infinite-dimensional counterparts but this requires verification 
using direct numerical simulations of the full PDE). Note that k > is the value of the 
imaginary part of the right-handed Bloch wall at its center [see Eq.flS])], and — k,q < is 
the imaginary part of the left-handed Bloch wall at its center. Therefore, the point (0, k,q) 
represents a configuration of the right-handed wall at x — +00 and the left-handed antiwall 
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at x = — oo. The point (0, — kq) corresponds to the left-handed wall at x — +00 and the 
right-handed antiwall at x — —00. 

Assume, first, that h is not close to or ~. When e is small and k is close to ±/«o (which 
are not small), terms in the third and second lines in ( 1361) are much smaller than the second 
term in the first line, and can be safely disregarded. The resulting dynamical system 

e = -^-A 2 €k(k 2 — Kq), k = 32heK 2 , 

does not have fixed points with nonzero e and k; on the other hand, the entire K-axis is a line 
of nonisolated fixed points. All these points describe pairs of walls moving with constant 
velocities ±i, where 

. _ 8(1 + h) 2 2 

Points on the positive-K semiaxis with k < k,q are stable and those with k > k,q unstable. 
Points on the negative-K semiaxis with k < — kq are stable and those with k > — k unstable. 
(See Figj9j) Consequently, walls with large initial separation [e(0) near 0] and i(0) > 
diverge to infinities (e — ► 0), whereas distant walls with i(0) < converge. 

As for the pairs of walls whose separation is not very large, the situation is straightforward 
for the walls which are not moving initially [i(0) = 0]; the corresponding initial conditions 
lie on the straight lines k = ±«o (dashed lines in Figj9j) These initial conditions give rise 
to the convergence of the walls (i.e. e grows as t 00). Otherwise, the type of interaction 
depends on whether the point (e(0), k(0)) lies to the left or to the right of the heteroclinic 
trajectory [the trajectory which connects the point (0, — k ) to (0, k ) — this trajectory is 
tangent to the vertical axis in Figj9j] 

These variational conclusions are in agreement with direct numerical simulations of the 
full PDE (J6j). (We employed the same numerical method as described in section IVD) The 
initial condition was chosen in the form of two Bloch walls of the opposite chirality which 
are initially at rest, i.e. Eq. lflO!) with K\ = K2 — Kq- The initially quiescent walls have always 
been observed to attract; see FigJlOl 

When h is close to |, k is small; since k is assumed to be in the vicinity of ±k , we 
cannot neglect terms in the third and second lines in ( 1361) in favour of the first line. However, 
it is sufficient to keep just the leading term (the term proportional to 7 — 9/i in the second 
line). A straightforward analysis reveals then that the situation here is similar to the one 
considered above: there are no fixed points with small nonzero e, while the K-axis is a line 
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of nonisolated fixed points, stable for k(k 2 — Kq) < and unstable otherwise. Therefore the 
phase portrait for small e and k is similar to the one above. 

When h is close to 0, the last term in the first line of (I36p becomes small and hence no 
terms in the third and second lines can be discarded. In this case, the dynamics can be 
analysed by plotting level curves of H using standard software. The same approach was 
adopted to study the vector field outside the small-e region. 

The resulting phase portrait is shown in Figj9]for h = 0.1; the portrait does not undergo 
any qualitative transformations as h is increased up to | or decreased down to 0. A striking 
feature of Figj9] is the existence of three families of closed orbits: one centred on a point 
with k ~ 1; a mirror-reflected family centred on a point with k « — 1; and a family of closed 
trajectories centred on a point on the horizontal axis. 

The closed trajectories centered on the fixed point with k ~ 1 represent a family of 
breather solutions of the NLS ([1]) (and so do their mirror- reflected twins). These breathers 
have been observed in our numerical simulations of the full partial differential equation 
(jEJ); see FigflQl They can also be constructed perturbatively (see section lIXl below). The 
nonlinear-center fixed point is sandwiched between two saddles, one below it (clearly visible 
in Figj9]) and one above. As h is decreased, the two saddle points approach each other. As 
a result, the periodic orbits are compressed in the vertical direction, and, as h reaches zero, 
the family of closed trajectories degenerates into a segment of the straight line k — 1. 

The closed loops about a nonlinear-centre fixed point with k = represent another family 
of temporally-periodic solutions. Since the value of k for a free-standing wall (i.e. k,q) is close 
to zero only when h is close to |, the family centred on n = admits a reliable interpretation 
only in this limit. When h = 0, the center point has the coordinate e ~ 0.06; as h grows 
to |, the centre moves towards the origin so that the family of closed orbits shrinks to 
a single point e = 0, n = 0. For h close to |, i.e. in the region where these solutions 
can be interpreted in terms of the two-wall Ansatz ( fTUI) . the closed trajectories describe a 
pair of Bloch walls of opposite chiralities executing periodic oscillations of their separation 
(but remaining far away from each other.) These oscillatory "doublets" have not yet been 
observed in direct numerical simulations and their existence remains an open problem. (See 
the concluding section for a possible explanation of the "nonvisibility" of these objects.) 
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VIII. TWO BLOCH WALLS OF LIKE CHIRALITIES 



Finally, we examine the case of a pair of Bloch walls (more precisely, a wall and an 
antiwall) of the same chirality. As in the previous section, we let 7 = and consider h < |. 
Assume, for definiteness, that we have a right-handed pair. When at rest, the well-separated 
wall and antiwall have equal (nonzero) momenta Pb = — vtk where k,q is as in (Q. When 
the walls start moving as a result of their interaction, the total momentum remains the same 
and hence the individual field momenta of the two walls will change by the same amount: 
one will become Pi = — tt(kq + q) and the other one P2 = — 7r(/? — q)- This means that 
one of the walls will have the amplitude of its imaginary part increase by q and the other 
one decrease by q, and so we will not have a symmetric configuration of an equal-amplitude 
wall and antiwall any longer. Therefore one may question the validity of our assumption 
that the two walls remain equal distance away from the origin for t > 0. To see that the 
equal-distance Ansatz fflOl) remains valid — at least for not too late times — we note that 
the velocity-momentum curve P(v) of a single Bloch wall has a nonzero slope at the point 
v = 0, P = Pb (see [5]; the curve is also reproduced in part II of the present series of papers 
[3]). This means that the wall whose momentum has increased by AP = nq, acquires the 
velocity Av, while the wall whose momentum has decreased by AP = irq, starts moving with 
the velocity —Av where Av = (dP/dv)~ 1 AP. Accordingly, if two walls are equal distance 
from the origin initially, they will remain equally far from the origin for all t > 0. This 
symmetric arrangement will break down only if AP is very large or if (dP/dv)^ 1 is very 
small, in which case one would have to take into account deviations of the shape of the P(v) 
curve from a straight line. The derivative (dP/dv)~ l tends to zero only if h — > | and indeed, 
numerical simulations do reveal nonsymmetric motion of like-chirality walls when h is very 
close to I (see below), but outside this narrow region the symmetric Ansatz (Tl0|) remains 
valid. 

Letting k\ = =F/to + q, ^2 = + Q an d B = B in the Lagrangian ( 1131) . we get equa- 
tions describing the dynamics of two like-chirality Bloch walls. These can be written as a 
hamiltonian system 

z= — q = (37a) 

TIT Oq TXT OZ 

where 

r(e, q) = 1 + 2e[l - 14e - 3(1 - 6e)«jj + (1 + 10e)g 2 ] > 



27 



and the Hamiltonian 



+ 3B* 



K 2 

7-9h- -£(7 + 38/i + 79h 2 ) 
A 2 



-4t(3 + I2h + h 2 )q 2 + (19 - h)q 4 - AA 2 q & 
A 2 

+ 16e 2 z[AhK 2 + 4(1 - h)q 2 - A 2 q% (37b) 

Assume, first, that h is not close to |, so that k is not a small quantity. In this case the 
Hamiltonian (137b 11 with small e and g is dominated just by two terms, 



H{z,q)*-±A*«tf-^«>e. (38) 

As in the previous section, the g-axis in this case is a line of nonisolated fixed points. 
Since z ~ — (16A 2 n 2 ] /3iiB)q, points with q > and q < represent pairs of converging 
and diverging walls, respectively. The trajectories are parabolas e = e* — (A 2 /6h)q 2 ; they 
start at points on the positive-g semiaxis and flow into points on the negative-g semiaxis. 
This implies that the initially diverging walls continue to diverge whereas the walls initially 
moving towards each other shall stop at some finite separation, and then repel to infinities. 
(A simple but practically important situation concerns walls which are initially at rest; these 
will obviously start repelling straight away.) The above conclusion applies, of course, only to 
slowly moving walls which are sufficiently far from each other [i.e. e(0) and g(0) are small]. 

Fig JTTT a) displays a typical result of direct numerical simulation of the interaction of two 
like-chirality Bloch walls within Eq. ([?)]). We tested several values of h and a variety of z(0). 
The initial condition was taken in the form (1TU1) with K\ = —k,q and k 2 = k,q. In agreement 
with conclusions of the variational analysis, the initially quiescent walls were observed to 
repel for all h not very close to |. (The case of h close to | is discussed below.) 

As h tends to zero, the phase portrait does not undergo any qualitative changes. However, 
the turning point of the parabola starting at a point g* on the g-axis, e* = (6/i) -1 / 2 Aq*, 
is shifted out of the small-e region, i.e. out of the range of applicability of variational 
approximation. Consequently, we can no longer expect the initially converging walls to stop 
a large distance away from each other and repel to infinities after that. Instead, for h ~ 
the colliding like-chirality Bloch walls should penetrate into the core of each other, and the 
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numerical simulations of the full PDE could be the only way to study the outcome of this 
deep impact. 

Next, let h be near ~. In this case the Hamiltonian (137bl) with small e and q reduces to 

H(z, q) « -^(4^ g 2 + g 4 + - 6eg 2 - 12e 2 ). 
9 

The dynamics with very small e and g is similar to the one described by Eq. (l38|) : distant 
pairs of walls with z > diverge while those which are initially set to converge (i < 0) slow 
down, stop and then repel to infinities [Fig Jl2T a)] . A new feature in the (h ~ |)-case is a 
saddle point at (e, q) = (jKq, 0) and a region of attraction which is separated from the region 
of repulsion by the stable manifold of the saddle [FigJT2Ta)]- This fixed point was previously 
discovered in a variational analysis of the (2:l)-resonantly forced Ginsburg-Landau equation 
[27I ] . Although it would be tempting to think that it represents a stationary complex of two 
Bloch walls of like chiralities, the actual interpretation of the stationary point turns out to 
be somewhat different. 
Indeed, it was proved in 



21] that in the absence of damping, two Bloch or two Neel 
walls cannot form stationary complexes. On the other hand, letting K\ = —k 2 = — an d 
e~ 2z = with small k , the trial function fflOl) can be written, approximately, as 



1m ijj 



e 2x + e 2x — e 2z 
e 2x + e~ 2x + e 2z ' 
A{e x - e~ x ) 

e 2x _|_ e -2x _|_ e 2z 



This coincides with the expression for the Bloch-Neel bound state {Eq.(ll) of Ref. 19J} with 



s = and a = —1, where we only need to send B — > 1 and make the identification e 2/3 = \e lz 



2 



Therefore, the "bound state" of two like-chirality Bloch walls produced by the variational 
analysis is, in fact, the Bloch-Neel complex with s = 0. The reason why the Bloch-Neel 
complex could be mistaken for a complex of two Bloch walls, was simply because the Bloch 
and Neel walls become indistinguishable as h — > |. 

The appearance of the region of attraction for h close to | deserves a special comment. 
The attraction becomes possible due to the smallness of k in this limit. Indeed, for suffi- 
ciently small kq, the momentum K\ = =p/«o + q becomes close to k 2 = ±Ko + q — which is 
characteristic for two Bloch walls of opposite chirality. Thus, in the limit h — > |, the walls 
may effectively change their chiralities. 
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It is instructive to consider the full phase portrait of the system (|37|) ; we produce it by 
plotting level curves of the Hamiltonian (137b j) beyond the region of small e and k [Fig JT^ b)]. 
A remarkable feature here is a family of closed orbits surrounding a nonlinear-center point. 
The center is born, along with the saddle point, in a saddle-center bifurcation as h is increased 
past 0.189. The presence of periodic orbits suggests that breathers may be formed in the 
collision of Bloch walls of like chirality. 

As h — > |, the saddle point 0) moves closer to the origin and so the periodic orbits 
pass near e = 0. Accordingly, in this limit the breathers should be accessible from initial 
conditions in the form of well-separated pairs of walls. We note, however, that all closed 
orbits pass also through a region of large e. [In fact, as h — > |, the coordinate of the 
nonlinear center tends to infinity: e c — > exp(l/KQ).] Consequently, the formation of the 
breathers requires accurate testing within the full PDE. In our numerical simulations with h 
close to |, we have indeed observed the attraction of walls with moderate separations (e.g. 
with e = 0.005 in the case of h = 0.33 and e = 0.001 in the case of h = 0.333) followed 
by the formation of a breather [Fig{TT](b)]. It is interesting to note that the attraction of 
two like-chirality walls with h ~ | proceeds via the chirality transmutation: first, one of 
the walls changes its chirality to the opposite; after that, the two opposite-chirality walls 
attract. 

Finally, we need to mention an anomaly in the interaction of two like-chirality walls which 
was observed in some of our simulations with h near |. The walls which are initially at rest 
and equal distance away from the origin, would start moving with different velocities or even 
in the same direction — violating our Ansatz (TTUT) based on the assumption of a symmetric 
arrangement for all t. This asymmetric anomaly is observed in a wider range of /i-values 
near h = | than the chirality transmutation (we detected the former for h as far from | as 
h = 0.30) and has a straightfoward explanation in terms of the P{v ) curve. As h — > ~, the 
derivative dP/dv tends to infinity at the points where the S"-shaped P(f)-curve intersects 
the P-axis and so the equation Av = (dP/dv)~ 1 AP has to be amended by keeping the 
symmetry-breaking term in (AP) 2 . It is this term that causes the asymmetric motion of the 
two walls. The anomalous interaction is typical for Bloch-Neel pairs and can be interpreted 
as the interaction of particles with the opposite mass sign ljj. 




The fact that the motion of two walls has to be asymmetric becomes even more obvious 
if one notices that the left and right turning points of the S*-curve approach v = as h — > ^. 
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As a result, the motion in the direction of the turning point becomes hampered. When h is 
extremely close to |, a new channel of interaction becomes available; namely, the wall which 
is repelled by its partner in the "hard" direction, may "tunnel" through the P(v )-barrier 
by means of the chirality transmutation, and the original symmetric arrangement remains 
undisturbed. 

IX. OTHER LOCALISED ATTRACTORS: NONDISSIPATIVE BREATHERS 

In the nondissipative situation (7 = 0), colliding walls form a moving or quiescent 
breather: a temporally oscillating, spatially localised structure. The breather was formed, 
for at least a short time, in every undamped collision that we simulated. The formation of 
the breather is accompanied by the release of a large amount of radiation; the numerical 
simulations reveal that the breather carries only about a quarter of the energy of the initial 
configuration, and so the remaining three quarters must escape into the radiation field. As a 
result, when the interval of simulation is too short, the interaction of the breather with the 
radiation re-entering the interval via the periodic boundaries is strong enough to destroy it 
after just a few oscillations. On a large interval, however, the amplitude of the re-entering 
radiation is small due to the dispersive broadening, and the breather persists indefinitely. 

The breather can be easily constructed perturbatively. First of all, a harmonic solution 
Sip = e 1 ^-^) of the NLS ([6]) linearised about ijj — 1 obeys the dispersion relation 

u 2 = M +2h + — J -1. 

The (k = 0)-harmonic oscillates with a nonzero frequency ujq = 2\/hA; therefore, one may 
also expect to find a broad small-amplitude breather oscillating with that frequency. We 
write 

if>(x,t) = l + e^ 1 (X 1 ,X 2 ,...,T ,T 1 ,...) 

+e 2 ^ 2 (A 1 ,A 2 ,...,T ,T 1 ,... ) + ..., (39) 

where X n = e n x, T n = e n t and e is a small parameter, and substitute (1391) in equation (jBj). 
Setting to zero coefficients of e 1 and e 2 , we find: 

r/>i = + 1P2 — U2 + W2, (40) 
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where 





l-l 


f 2h\ 









^-e M + c.c. 
2h 



and 





l-l 


' l + 2/i\ 






1 iwo / 



2iwoTo 



2/i 




+ c.c. 



(41) 



(42) 



Here the amplitude depends only on long scales X\, X 2 , ... and slow times 7\, T 2 , while 
c.c. indicates the complex conjugate. The solvability condition at the order e 2 is d<p/dTi = 0, 
so that 0, in fact, is independent of T\\ = 4>(X 1 , X 2 , . . . , T 2 , T 3 , . . . ). 

At the order e 3 , the solvability condition forces to obey the attractive unperturbed 
nonlinear Schrodinger equation 



-Ai 



o. 



(43) 



The breather results if we choose in the form of the soliton 

= 0o sech 



1+ Ah 



x exp 



l + 2h 



A(l + Ah) Vh 



(2X, - VT 2 ) 



(44) 



where 0o and V are free parameters describing the amplitude and velocity of the soliton, 
respectively. Substituting Eq.( l44l into ( !39l) - (l4B and defining the amplitude of the breather 
a = e0o and its velocity v = eV, we obtain, finally, 



Re ip = 1 + 2a sech 



Im-0 = —sech 



. 1 + Ah , . 

2\ -a(x — vt) 

1 + 2h v ' 



Vh 



,1 + Ah. . 
2\l —a(x — vt) 



l + 2h 



cos 



sin 



Qt-v 



Qt-v 



Vh 



A(l + 2h) 



(2x - vt) 



Vh 



A(l + 2h) 



(2x - vt) 



+ G(e 2 ),(45a) 



+ 0(e 2 ), (45b) 



where 



A{l + Ah) 2 

Q = uoq — — — -a 1 . 

yh 



(46) 



We carried out numerical simulations of Eq.Q with Eq. (1451) as an initial condition, with 
a variety of (small) values of a. (We confined our simulations to the case v — 0.) For 
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all values of h and a we tried, the breather persisted for the full length of the simulation 
(approximately 1000 periods of the breather's oscillation), with virtually no change. The 
measured frequency of oscillation coincided with the asymptotic value (jlBjl up to the third 
decimal place. 

Simulations were also performed with a gaussian initial condition 

^ = l + ae- px \ (47) 

for a number of small complex values of a (\a\ < 0.5) and positive f3. All runs resulted in 
the formation of a breather, although in some cases [for \a\ > 0.3], the emerging breather 
would break into a pair of counterpropagating breathers. In those cases where the emerging 
breather was nonmoving, its shape and frequency were found to be in an excellent agreement 
with the asymptotic formula (H5]) . 

X. CONCLUSIONS AND OPEN PROBLEMS 

In this paper, we studied the interactions between the similar-type dark solitons of the 
nonlinear Schrodinger equation, i.e. Bloch-Bloch and Neel-Neel interactions. Our approach 
was based on the variational approximation (also known as the collective-coordinate, or 
variation-of-action, method). The variational conclusions were verified using numerical con- 
tinuation of solutions of the stationary damped-driven NLS fl33l and via direct numerical 
simulations of the full partial differential equation ([6]). 

A. Conclusions 

In the dissipative situation (7 ^ 0), the only available solitons are the Neel walls. When 
two Neel walls are very far apart, their interaction is simple: the walls repel if h 2 < 7 2 + | 
and attract if h 2 > •y 2 + |. The repelling walls diverge to infinities; as for the case of 
attraction, there are two possible scenarios. In order to distinguish between the two, one is 
led to consider the situation where the walls are closer to each other [but still sufficiently 
far apart for the large-separation approximation (fit)]) to remain valid.] 

At these shorter distances, the dynamics is influenced by two bound states, a stable and 
an unstable one. The stable bound state exists for h between ^Jy 2 + 1/9 and a threshold 
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driving strength h sn , and the unstable one exists for all h < h sn . Here h 2 sn = 7 2 + (A 2 sn — l) 2 
where A 2 sn {^f) is defined as a root of the system (I3~i~l) - fl32l) . For small 7, the curve h sn (j) can 
be described explicitly: 

h sn = i + 0.3048 7 2/3 . (48) 

[This is Eq. (l29[) written in terms of h and 7.] In their region of coexistence, the stable 
complex has a larger separation: 2z s > 2z u . 

When h is smaller than ^7 2 + 1/9, the walls repel if their separation distance 2z(0) is 
greater than 2z u , the interwall separation in the unstable complex. If 2z(0) < 2z u , the 
walls converge and annihilate. (This verdict does not extend to the region where both h 
and 7 are small. In this region the variational analysis of the small-separation dynamics is 
inconclusive.) When ^7 2 + 1/9 < h < h sn , pairs of Neel walls with separations 2z(0) larger 
than 2z u evolve towards the stable bound state while those with 2z(0) < 2z u converge and 
annihilate. Finally, the walls with h > h sn converge and annihilate irrespectively of their 
initial separation. These results pertain to walls at shorter distances (which are however 
sufficiently far apart for the variational approximation to remain valid). In particular, they 
answer the question as to what finally happens to the two walls attracted from very large 
distances. 

The nondissipative case requires a separate summary. Here the walls can move at constant 
speeds and the interaction pattern becomes complicated by the presence of inertia. When h 
is greater than | (by small or large value), the Neel walls attract and converge — unless the 
initial condition corresponds to walls having large and opposite velocities. In the latter case 
the attraction is unable to stop the diverging walls and they escape to infinities. On the 
other hand, when h is smaller than |, the walls repel. The exception here is the case where 
h is close to |; in this case walls with very large separations repel whereas walls which are 
not so far from each other, attract. 

In the dissipation-free case, the available dark solitons also include Bloch walls. The 
interaction between two Bloch walls depends on their relative chiralities: two initially qui- 
escent, oppositely-handed Bloch walls attract while two quiescent walls with like chiralities 
placed at a large distance away from each other, repel. The exception is the case of h close 
to |; in this limit, two walls of like chirality repel at large distances but exhibit anomalous 
interaction or transmute into an opposite-chirality pair and attract - when placed closer to 
each other. 
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These conclusions can be extended to the case of the moving Bloch walls, where one just 
needs to take their inertia into account. For example, two initially diverging oppositely- 
handed walls at large separation will continue to diverge despite the attraction whereas two 
likely-handed walls which were initially moving against each other, will continue to converge 
(until the repulsion stops them and sends away to infinities). 

In addition to the interactions between well-separated walls, we investigated products 
of their collision. When the system is not damped, the collision of two walls results in a 
stationary or travelling breather. We reconstructed the numerically found breather as an 
asymptotic series. 

When 7 is nonzero, oscillations are damped and the only nontrivial product of collision 
of two Neel walls is a stationary bubble — the bound state of two Neel walls. Using the 
numerical continuation of solutions to the ODE fl33l) . we demarcated the bubble's domain 
of existence in the parameter space. [This domain is, naturally, a subset of the part of 
the (7, /i)-plane in which remote walls attract]. The numerical demarcation is in excellent 
agreement with the domain of existence obtained variationally. 

B. Open Problems 

1. We complete the paper by listing several open problems to which we plan to return in 
future. One interesting problem that merits further investigation, concerns the (unstable) 
bubble solution with small h and 7. (This solution is exemplified by FigfTBl) Our numerics 
and the variational analysis in section IIVDI show that this solution cannot be regarded as 
a complex of two Neel walls, not even strongly overlapping ones. Its asymptotic behaviour 
is closer to that of a Bloch wall; however, Bloch walls do not exist in presence of damping. 
Furthermore, the bubble cannot be continued to 7 = and hence does not have a nondis- 
sipative analogue in the form of a bound state of two Bloch, or a Bloch and a Neel, walls. 
(This follows from the fact that the Neel wall is the only solution in the 7 = domain which 
admits continuation to nonzero 7 as it is the only solution with zero momentum [5j.) Thus a 
question arises whether the bubble with small h and 7 could be interpreted as a complex of 
two hypothetical "dissipative Bloch walls" of opposite chiralities — which do not exist indi- 
vidually but can exist as a bound state due to the cancellation of their opposite momenta. A 
similar momentum cancellation occurs in the attractive damped-driven NLS where it allows 
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complexes to travel with nonzero speeds despite the fact that their constituent solitons are 



immobile 
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A related question concerns the short-distance Neel wall dynamics with small h and 7. 
This problem is inaccessible to the variational method, and the direct numerical simulations 
of the full PDE seems to be the only appropriate line of attack here. 

2. With regard to the Bloch walls, an open problem concerns the family of closed orbits 
centred on tz = in Figj9l When h is close to |, these periodic orbits can be interpreted 
within the full partial differential equation (pQ), as pairs of opposite-chirality Bloch walls 
executing periodic oscillations of their separation. It is not impossible that these periodic 
solutions exist also for h not so close to ~ [where the validity of the well-separated Ansatz 
( FlOl) becomes questionable.] It remains to be understood why these orbits have not been 
observed in our numerical simulations of Eq.([T]), with either h. One possible explanation 
stems from the fact that the nonlinear-center fixed point enclosed by this family of closed 
orbits, corresponds to a maximum of the energy ( 113cl) . Consequently, the nonlinear modes 
not captured by our two-dimensional Ansatz (flOl) — such as radiations — should make 
the point (e, k) slide downhill towards orbits with larger radii, until it crosses through the 
homoclinic trajectory and ends up at e = 0. [It is pertinent to note here that the centre 
points enclosed by the other two families of closed orbits in Figj9] are minima of the energy 
(I13cp ; this is consistent with the stability of the corresponding breather solutions observed 
in numerical simulations in Section IVTl ] 

3. Next, it would be interesting to simulate collisions of the like-chirality Bloch walls 
with h ~ 0. This case is not amenable to the variational analysis due to the large extent 
("width") of the walls. 

4. It would also be interesting to gain a deeper theoretical insight into the anomalous 
interaction of the like-chirality Bloch walls arising for h close to ~. The variational ap- 
proach remains a proper tool here; one should only generalise the Ansatz ( flOl) by allowing 
nonsymmetric configurations of the walls: 

<fi(x, t) = tanh [B(x + Zx)) — iKisech[B(x + z\)], 
f 2 (x,t) = tanh [B(x — z 2 )} + m 2 sech[i3(x — z 2 )}- 

Here zi, z 2 and K\, k 2 are unrelated pairs of variables. 

5. Finally, one more future challenge is the analysis of the interaction of two breathers 
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and their synchronisation to a common frequency. 
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FIG. 3: The vector field (125!) . (a) fj, < (small initial conditions are attracted to a stable fixed point at 
the origin; larger initial conditions give rise to trajectories flowing to infinity), (b) < fj, < 0.30487 2 / 3 
(initial conditions with small e and k are attracted to a stable fixed point away from the origin; larger 
initial conditions give rise to unbounded trajectories), (c) /i > 0.30487 2 / 3 (no stable fixed points, all 
trajectories flow to infinity). The stable fixed point is marked by a solid blob; the unstable one by an 
open circle. 
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FIG. 5: A bubble with a small distance between the walls (a tightly bound complex) looks like a single 
entity. (Note that the real part does not even reach —1 in the core region.) The solid and the dashed 
lines stand for the real and imaginary part of if), respectively. For comparison, we also show a bubble 
with a wide interwall separation arising for h 2 close to 7 2 + l/9 (inset). Here only the real part is shown 

as Imtp(x) = 0. 
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FIG. 6: The energy of the dissipative bubble for the fixed 7 (main panel) and fixed h (inset). The solid 
and dashed branches indicate stable and unstable solutions, respectively. The stable branch terminates 
at the point where the distance between the walls becomes infinite; this point is found to coincide with 
A 2 = |. This termination scenario was predicted in section [TVl using the variational method. 
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FIG. 7: The existence and stability diagram of the bound states on the (7, /i)-plane (obtained numeri- 
cally.) The shaded corridor marks the region of coexistence of the stable and unstable bubble. Its upper 
boundary corresponds to the turning point where the two bubbles merge; there are no bound states 
above this line. The dotted curve gives the variational approximation ()31[) - t|32[) to the line of turning 
points. The lower boundary of the shaded area corresponds to the stable bound state of two infinitely 
separated walls (an infinitely long stable bubble). This numerical curve is visually indistinguishable from 
h = \/7 2 + 1/9. In the strip between h = yT + 1/9 and h = 7, only the unstable bubble is found. 
Inset: a blow-up of the small-7 region. Here, the dotted curve is plotted using an asymptotic formula 
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FIG. 8: Attraction of two dissipative Neel walls may result either in the annihilation of two walls (a) 
(here 7 = 0.1, h = 0.5) or formation of a bound state (c) (here 7 = 0.1, h = 0.41). In the absence of 
damping, the attraction results in the formation of a stationary breather (b) (here h = 0.5). Simulations 
reported in panels (a,b) correspond to a point (7, h) lying above the shaded region in Fig|7]while those 
in panel (c) pertain to a point inside this region. All panels show the real part of ift, multiplied by — 1 
for better visualization. 
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FIG. 9: Trajectories of the system (}35j) obtained as level curves of the Hamiltonian (1361) . (Here 
h = 0.1.) Two points at which the trajectory is tangent to the K-axis, correspond to stationary pairs 
of walls. The dashed lines are k = ±ko; the initial conditions lying on these lines pertain to stationary 
walls: i(0) = 0. 




FIG. 10: Attraction of Bloch walls with opposite chirality; this attraction results in a long-lived breather. 
(In this plot, h = 0.1). To highlight the chiralities of the Bloch walls, the imaginary part of i/j is shown. 
(We remind that in the case of the wall-antiwall configuration that we are considering, the same signs 
of the imaginary parts imply opposite chiralities of the walls.) 
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FIG. 11: Interaction of Bloch walls of like chiralities. The imaginary part of t/j is shown, (a) Repulsion 
for generic h. (Here h = 0.15). (b) For h close to ^, the walls may attract and produce a breather on 
collision. (In this plot, h = 0.333). Note that before the walls start converging, the wall on the left 
changes its chirality to the opposite. (We remind that in the case of the wall-antiwall configuration 
that we are considering, opposite signs of imaginary parts imply like chiralities of the walls.) 




FIG. 12: Trajectories of two Bloch walls with like chiralities on the (e,g)-plane. (a) h close to | 
(h = 0.33 in this plot). Note that the saddle point has small e here; as a result, there is a region of 
attraction within the range of validity of the variational approximation, (b) Generic h (h = 0.2 in this 
plot). The family of closed orbits describe breather solutions of the NLS ([1]). Note that unlike in (a), 
the saddle point and the region of attraction arise at small interwall separations in (b). 
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FIG. 13: The tightly bound complex of two Neel walls with small h and 7. The solid and dashed line 
depict the real and imaginary part, respectively. Note an enormous spatial extent of the bubble. 
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